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Abstract. We study inflationary universe models that are characterized by a single scalar 
inflaton field. The study of these models are based on two dynamical equations; one corre- 
sponding to the Klein-Gordon equation for the inflaton field, and the other, to a generalized 
Friedmann equation. After describing the kinematics and dynamics of the models, we de- 
termine in some detail scalar density perturbations and relic gravitational waves. We apply 
this approach to the Friedmann-Chern-Simons and the brane-world inflationary models. 
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1 Introduction 

The early accelerated expansion of the universe (inflation) has been proposed as a good ap- 
proach for solving most of the cosmological "puzzles" [1, 2]. This brief accelerated expansion 
serves, apart of solving most of the cosmological problems, to produce the seeds for the large 
scale structure formation. In fact, the present popularity of inflation is entirely due to its 
ability to generate a spectrum of density perturbations which lead to structure formation in 
the universe. In essence, the conclusion that all the observations of microwave background 
anisotropy performed so far support inflation, rests on the consistency of the anisotropy with 
an almost Harrison-Zel'dovich power spectrum predicted by most of the inflationary universe 
scenarios [3]. 

The implementation of inflationary models rests on the introduction of a scalar inflaton 
fleld, (j). The evolution of this field becomes governed by its scalar potential, V{(j)), via the 
Klein-Gordon equation. In this way, this equation of motion, together with the Friedmann 
equation, form the most simple set of field equations, which can be used to study inflationary 
solutions. 

The previous implementation works if we give an explicit expression for the scalar infla- 
ton potential, V^(^). However, in most cases result complicated to flnd solutions, even in the 
situation in which it is applied the so-called slow-roll approximation, where the kinetics terms 
is much smaller than the potential energy, i.e. <C V{(p), together with the approximation 



Another way to find inflationary solutions, out of the slow-roll approximation, is giving 
the functional form of the Hubble parameter in term of the inflaton fleld, i.e. -ff(</>) [4]. At 
first glance this approach looks quite unsuitable, since the function -ff(</)) could be chosen 
arbitrarily. Of course not any Hubble function, H{(j)), will yield an appropriated infiationary 
solution. The same it is found when the slow-roll approximation is taken into account. 
There, specific scalar potentials are chosen in order that the model could be implemented. 
For instance, if the scalar potential does not present a minimum, then, the exit from infiation 
becomes a problem, since the usual reheating process can not be carried out. 

It seems that the description of inflation in term of the Hubble parameter as a function 
of the scalar fleld, i.e. H{(j)), looks "more natural" [5]. In this way, this approach known as 
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the Hamilton- Jacobi formalism (H-J) provides us with a straight-forward way of exploring 
the inflationary scenario [6, 7]. 

The H-J approach presents some advantages when compared with the slow-roll approx- 
imation: first, in the H-J formalism the form of the potential is deduced, and second, since 
an exact solution is obtained, then, application to the final period of inflation is possible, 
where the kinetic term of the inflaton field in the Friedmann equation becomes important, 
i.e. when studying the final stage of inflation. In this approach we can use the scalar field, 
(j), as a "time variable", and for that, we demand that this field increases monotonically, i.e. 
its time derivative, should not change of sign along the inflationary evolution. 

In this article we would like to study the consequences that result when considering a 
modified Friedmann equation expressed by 
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(1.1) 



where > is an arbitrary function of the Hubble parameter H = a/ a, with a the scale 
factor, the prime represents a derivative with respect to the scalar inflaton field 4), i.e. ' = ^ 
and the dots represent derivatives with respect to the cosmological time, t. Here, the inflaton 
potential is expressed by V{4)) and m?pi = 1/G represents the Planck mass. 

The motivation for using this kind of equation lies in the fact that in the literature 
several models have been studied which can be reduced to such modified Friedman equation. 

2 

In a L(i?)-theory of gravity in which L{R) = R — where R is the scalar curvature and a 
is a constant with dimension of mass square, the Friedmann equation becomes modified by 

the expression J^{H) = -r^ g — ^[8]- 



Also, it is possible to consider ^{H) = H — aH , where a is a constant with di- 
mension of [mass]~'^. There exist various forms in arriving to this expression for J^{H). 
This has been derived by considering a quantum corrected entropy-area relation of the type 

Sa = ^'pi'^ — a In (y^'pi^^ t^]' where A is the area of the apparent horizon, and a is a 

dimensionless positive constant determined by the conformal anomaly of the fields. This 
conformal anomaly is interpreted as a quantum correction to the entropy of the appar- 
ent horizon[10]. On the other hand, this modified Friedmann equation could be obtained 
when an AdS-Schwarzschild black-hole via holographic renormalization is considered, to- 
gether with mixed boundary conditions corresponding to the Einstein field equations in four 
dimension[ll]. Also, this could be derived in terms of spacetime thermodynamics together 
with a generalized uncertainly principle of quantum gravity [12]. A Chern-Simons type of the- 
ory yields to this modification too[13]. The resulting Friedmann equation when considering 
this type of modification we will call the Friedmann-Chern-Simons equation[14]. 

On the other hand, superstring and M-Theory bring the possibility of considering our 
universe as a domain wall embedded in a higher dimensional space. In this scenario the 
standard model of particle is confined to the brane, while gravitation propagate into the 
bulk spacetime. The effect of extra dimensions induces a change in the Friedmann equation. 

SttA 



Here, the function J-{H) results to be J'{H) 



1 + 
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where A 



represents the brane tension[15]. 

Here, in this paper, after giving a general approach to the study of infiation based on 
the modified Friedmann equation, we will describe in some detail the latter two cases speci- 
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fied previously, i.e. the Friedmann-Chern-Simons and the brane-world inflationary universe 
models. 



2 The exact solution approach 

Our study will be based on considering Eq. (1.1) together with the scalar field equation 

4> + 3H,i) + V'{(p) = 0. (2.1) 

In obtaining this latter equation we have assumed that the matter, specified by the inflaton 
scalar field, enters into the action Lagrangian in such a way that its variation in a Friedmann- 
Robertson- Walker background metric leads to the Klein-Gordon equation, expressed by Eq. 
(2.1). Therefore, we are considering constrained sort of models, in which the background 
(together with the perturbed equation, see Eq. (3.3)) are not modified. In this context 
theory of gravity , such that IIofava-Lifshitz[16] lies outside of the approach followed here. 
Also, in this study we will use the scalar field i;^ as a "time variable" . The requirement imposed 
in this approach is that the scalar field increases monotonically and its time derivative, (j), 
should not change of sign along the path evolution. 

From the field Eqs. (1.1) and (2.1) it is found that 
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where J-"^ = dJ- /dH. This latter equation allows us to write down an explicit expression for 
the scalar potential 



y(0) = 

It is not hard to show that 
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from which we get 
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where = a{(j).). 

On the other hand, we see that the acceleration quation for the scale factor results to 

be 
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where the function corresponds to 
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From this latter expression we can see that this definition, called the first Hubble slow-roll 
parameter, gives information about the acceleration of the universe. During inflation we have 
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that < 1, and this period ends when takes the value equal to one. In the next section 
we will use this parameter for describing scalar and tensor perturbations. 

One interesting quantity in characterizing inflationary universe models is the amount of 
inflation. Usually, this quantity is deflned by 

a{te) 
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where a (te) corresponds to the scale factor evaluated at the end of inflation. 
Friedmann Equation it becomes, in terms of the scalar field 
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For the modified 

(2.9) 



Here, (pe represent the value of the scalar field at the end of infiation. Its value is determined 
by imposing that {(pe) = 1- 

It seems to be more appropriated to describe the amount of inflation in terms of the 
comoving Hubble length, l/{aH) than in terms of the scale factor only. In this case the 
amount of inflation becomes deflned as[17] 

a{te) H{Q 
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which results into 
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Note that, in general, N((/)) is smaller that A^(i;^) and only in the slow-roll limit they coincide. 

In the description of inflation it is convenient to show that their solutions are indepen- 
dent from their initial conditions. This ensure the true predictive power that presents any 
inflationary universe model, otherwise the corresponding physical quantities associated with 
the inflationary phase, such that the scalar or tensor spectra, would depend on these initial 
conditions. Thus, with the purpose of being predictive, any inflationary model needs that 
their solutions present an attractor behavior, in the sense that solutions with different initial 
conditions should tend to a unique solution [5]. 

In order to study the corresponding inflationary attractor solutions for our case, we 
follow Ref. [5]. We start by considering a linear perturbation, 6H{(j)), around a given infla- 
tionary solution, expressed by i?o (</>)• In the following we will refer to this quantity as the 
background solution, and any quantity with the subscript zero is assumed to be evaluated 
taking into account the background solution. Therefore, at flrst order on 5H{(j)), we get from 
the field Equations (1.1) and (2.1) that 
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This latter expression can be solved for getting 
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where (p. corresponds to some arbitrary initial value of (j). By considering theories in which 
J-"^ > H T we find that the integrand within the exponential term will be negative, since 
and H' have opposite signs (assuming that cj) does not change sign due to the perturbation 
5H). Thus, all the linear perturbations tend to vanish quickly [17]. 
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3 Scalar and tensor perturbations 



Inflation generates perturbations through the amplification of quantum fluctuations, which 
are stretched to astrophysical scales by the brief, but rapid inflationary expansion. The 
simplest models of inflation generate two types of perturbations, density perturbations which 
come from quantum fluctuations of the scalar fleld[18-21], and relic gravitational waves which 
are tensor metric fluctuations [22-26]. The former experience gravitational instability and 
lead to structure formation[27], while the latter predicts a stochastic background of relic 
gravitational waves which could influence the cosmic microwave background anisotropy via 
the presence of polarization in this [28]. The upcoming experiments such as Planck satellite 
will characterize polarization anisotropy to a higher accuracy [29]. It is very timely to develop 
the tools which can optimally utilize the polarization information to constrain models of the 
early universe. Specifically, the magnetic modes (B-modes) are signal from cosmic infiation 
and suggest the presence of gravitational waves [30]. 

In order to describe these perturbations let us consider the Hamilton- J acohi slow-roll 
parameters. The first slow-roll parameter was already defined in the previous section, Eq. 
(2.7). 

The second slow-roll parameter, r/^, is defined as 

dlna ~\87r J H H- ^"^'^^ 

We should note here that both and r]^ are exact quantities, although we call them "slow- 
roll parameters". In the slow-roll limit these parameters becomes[17, 31] 

Vh — > V-^, (3-2) 

where the quantities e and r] are the common slow-roll parameters which satisfies e ^ 1 and 
r] <^ 1, in agreement with the slow-roll approximation. 

The evolution equation for the Fourier modes of the scalar perturbations (quantum 
mode functions) at some comovil wave number scale k is governed by[32, 33] 

-^+[^ --^J-^ = 0^ (3-3) 

where r represents the conformal time defined by r = j \dt and corresponds to the 

Fourier transformed of the Mukanov variable, which is defined hy u = zTZ, with z = ajj and 
TZ defining the gauge-invariant comovil curvature perturbation. This latter amount remains 
constant outside the horizon, i.e. metric perturbations with wavelengths larger than the 
Hubble radius will be frozen[34]. 

During infiation it is expected that fe^ ^ physical modes are assumed to 

have a wavelength much smaller than the curvature scale, and thus Eq. (3.3) can be solved 
to achieve 

nfc(r)~e-'^-fl + ^ + ....V (3.4) 



On the other hand, when hP' <C ^^7!) we have that the physical modes correspond to wave- 
lengths much bigger than the curvature scale. 
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The mass term -4-4 becomes in our case 



For = H"^ we obtain that [35] 
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It has long been known that Eq. (3.3) could be solved exactly in the case in which 
the mass term is proportional to r~^, in which case this equation reduces to a Bessel 
equation, where the standard solution becomes Uk ~ \/ —kTH^{—kT), with if^ the Hankel 
function of first kind, and the parameter u depends on the slow-roll parameter e via v = 
3/2 + e/(l — e). For instance, this occurred in the case of standard Friedmann equation and 
the scale factor, a(t), expands as a power law, i.e. a{t) ^ {p > 1), in which case it is 
obtained that ~ Vh ~ Constant [36]. Others solutions which are far from the slow-roll 
approximation are described in Ref. [37]. 

Immediately we get an explicit expression for u^, we can obtain the power spectrum, 
which is defined in terms of the two point correlation function as 

Vnik) = ^< '^t''^t > + (3.6) 



which in terms of the Uk and z it becomes 
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Vm = ^_\^\- (3.7) 

In order to obtain Uk by solving equation (3.3), we need to impose some boundary 
conditions. These asymptotic conditions are usually taken to be the so-called Bunch-Davies 
vacuum states [38] 

{-^e~'^^'^ as -kri — > 00, 
t ^ n (^-^^ 
AkZ as —KT] — > U. 

This ensures that perturbations that are generated well inside the horizon, i.e. in the region 
where k <C aH, the modes approach plane waves and those that are generated well outside 
the horizon, i.e. in the region where k ^ aH, remain unchanged. 

The description of the primordial curvature perturbation presents a standard result 
given by [20, 21, 
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(3.9) 

aH=k 

This perturbation is in general a function of the wave number k, which is evaluated for 
aH = k, i.e. when a given mode crosses outside the horizon during inflation. Since the 
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modes do not evolve outside the horizon, the amphtude of the modes when they cross back 
inside the horizon, coincides with the value that they had when they left the horizon. 

By using the primordial scalar perturbations we can introduce the scalar spectral index 
Us defined by 

_ dhiVn 

This quantity becomes 

ns-l = 2r^^-2[?,-H^^e^. (3.11) 

In the same way we define the running scalar spectral index, ag = — — which results to be 
given by 

a, = 2 ("s - m^^\ e^ri^ - 2 U - + H^^^] el -2(^, (3.12) 

where corresponds to the third slow-roll parameter and becomes given by 

( ^ ) (3.13) 

In addition to the scalar curvature perturbations, transverse-traceless tensor pertur- 
bations can also be generated from quantum fluctuations during inflation [23, 27]. The 
tensor perturbations do not couple to matter and consequently they are only determined by 
the dynamics of the background metric, so the standard results for the evolution of tensor 
perturbations of the metric remains valid. The two independent polarizations evolve like 
minimally coupled massless fields with spectrum^ 
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(3.14) 

aH=k 



Similarly to the case of scalar perturbations, we evaluate the expression on the right hand 

side of Eq. (3.14) when the comoving scale k leaves the horizon during inflation. Furthermore, 

dlnVr 

dink 
results to be 

n^ = -2e„. (3.15) 



we can introduce the gravitational wave spectral index n^, defined by n^, = — - — —, which 



Here, we can also introduce the running tensor spectral index defined by 

At this point we can define the tensor-to-scalar amplitude ratio r = ^— which becomes 

^ = 2^e„ (3.17) 



^We mention here that this expression should be implemented with a factor such that F^{H/fi), where 
Fa^i^) = Vl + ~ ( i+4ap^ ) •''^ sinh~^ i, when a braneworld with a Gauss-Bonnet term is considered [40]. 
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and combining equations (3.15) and (3.17) we find that 




(3.18) 



This latter expression corresponds to the inflationary consistency condition[35 , 41]. Note that 
for standard cosmology, in which T{H) = H^, Eq. (3.18) reduces to r = —i^n^. However, 
this relation could be violated in some cases [42, 43]. Note that this relation depends on the 
kind of theory that we are dealing with. 



4 The hierarchy of the slow-roll parameters and the flow equations 

There exist a different way of studying inflationary universe models, which is subtended by a 
sort of hierarchy imposed on the slow-roll parameters [44, 45]. In fact, the set of equations in 
this approach is based on derivatives with respect to the e-folding number over the slow-roll 
parameters. 

We have previously introduced the slow-roll parameters, such as e^j , ry^^ and , to which 
we have given a sort of hierarchy calling them first, second and third slow-roll parameters, 
respectively. Each of these parameters is characterized by their dependence on the order 
of the scalar field derivative of the Hubble ratio, H{(p), such as ~ {H')'^, rj^j ~ H" and 
~ H'" , as we can see from equations (2.7), (3.1) and (3.13), respectively. It is possible 
extend this definition to higher derivatives of the Hubble parameter so that we can introduce 
the following parameter 



m 



in 



H 

w 



(4.1) 



where for Z = 1 we have that ^X^^ = rjj^ and I = 2 corresponds to ^A^ = 
It not hard to show that the following set of equations is satisfied 
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Here, it was used the relationship 
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In the standard case, i.e. when F{H) = H^, the above set of equations reduces 
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where a = 2rijj — Ae^j [45] . 
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In order to solve the infinite set of equations (4.2) the series is truncated by imposing 
a vanishing value to a given high enough slow-roll parameter. This corresponds to take 
that ^''"*'^A^ = 0, for an appropriated large number M (for instance, in the literature has 
been used M = 5 [45]). With this truncation the set of equations has been solved both 
numerically [44-47] and analytically [48-50]. 

With respect to possible solutions of these equations and their relations with the infla- 
tionary paradigm, it was emphasized in Ref. [48] that there is no a clear connection between 
them. Actually, it is not clear that a particular solution of the flow equations corresponds 
directly to some type of inflationary solution. To achieve this task we must add additional 
ingredients to the corresponding solutions. The main ingredient that has been left out of this 
scheme has been the Friedmann equation itself. Thus, in solving the flow equations we could 
get as a function of the scalar field (j) (imposing the condition that this parameter will 
satisfy the range ^ ^ 1), and then we could get H{(f)) through the following relation 



Thus, in principle, we could give the function F{H) so that we could get as a function of 
(p by solving the set of Eqs. (4.2) we could obtain the Hubble parameter, H((p), as a function 
of the scalar field through Eq. (4.4)^. With this in hand, we can obtain an explicit expression 
for the scalar potential, V{cj)), given by 
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In short, for a given function T{H), we could say that the Hubble flow formalism allows 
us to determinate the scalar field potential, V{(j)), associated to some inflationary universe 
model. In order to realize this task we first solve the flow equations, Eqs. (4.2), from which 
we can obtain the first slow-roll parameter, , under the condition that this parameter must 
satisfy the bound ^ ^ 1 . Then, by using Eq. (4.4), we get the corresponding Hubble 
parameter as a function of the scalar field, -ff(0), from which we could obtain all the other 
quantities associated to the inflationary scenario. 

5 Two interesting cases 

5.1 The Priedmann-Chern-Simons model 

As stated in the introduction we would like to consider here a model in which the Friedmann 
equation modifies to"^ 

rm^H^-aH^^iJ^y,, (5.1) 



^In the standard case in which = it is obtained that 

H[cj>) = exp 



Stt 



^For more details on this case see Ref. [14]. 
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where a is an arbitrary constant with dimension of . Here, we assume that during the 
inflationary evolution the Hubble parameter, i7, satisfies the bound K < 1/y/a, so that the 
energy density associated to the scalar field (p is positive. 

It is possible to choose for the generating function a polynomial like H{(j)) = Hq{1 + 
+ P24''^ + ■■■■(^n4'^), where Hq and the different /? are constants. This sort of solution was 
used to generate suitable functions of slow-roll parameters [48]. Here, just for simplicity, and 
in order to show how the this approach work, we shall take the previous polynomial, but, up 
to first order in the scalar field (p, i.e. H((j)) = Hq{1 + /3(p), with /? an arbitrary constant with 
In this case, the scalar potential becomes 
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where (/> = 1 + /3</). 

In the slow-roll approximation, i.e. where 
it is found that the scalar potential becomes 
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Figure 1 depicts the shape of the potential for the exact case (thick line), expressed by Eq. 
(5.2), together with the approximated slow-roll case, Eq. (5.3). In the same figure the dotted 
line represents the exact case in which the a-parameter is vanished. 
The scalar field results to be given by 
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where Ti = V2a {l3 Hq) 
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and (/'(to) = 0. This latter expression allows us to write down 



the Hubble parameter as a function of time. From this result we get the scale factor, a{t), 
which results to be 



a{t) = ao 



^sinh {2 tanh"^ [tanh (i cosh" W2a i^o)) e 



sinh [cosh ^{\/2a Hq 



(5.5) 



In order to see if this latter expression describes an accelerated phase, for given values 
of the parameters, we plot in Figure 2 the deceleration parameters q, which is defined as 
? = ■ For this plotting, we have taken the value ^^/aHo = -g^- The different curves 



correspond to different values of the exponent that appears in the scale factor a, i.e. 4 v 2 a 
These curves show that the universe is accelerating, since the parameter q(t) turns out to be 
negative as time passes. Therefore, our model presents a period of inflation, at least for the 
values of the parameters that we have considered here. 
The amount of inflation becomes in this case 
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Figure 1. Plots of the scalar potentials, V{^), as a function of the dimcnsionless "scalar field", $ = 
Hq (f>. The thick line represents the exact potential, expressed by Eq. (5.2). Dashed line represents 
the same potential, but in the slow-roll approximation, Eq. (5.3). The dotted line corresponds to 
the exact case, but when a = 0. Here we have taken a (/?i?o)^ = g^^? and V{^) is expressed as a 

multiple of the constant Vn = J"-^' . 




X 

Figure 2. Time evolution of the deceleration parameter, q, as a function of x = "H (t — to). Here we 
have taken the value y/aHt^ — -g^- The thick, dotted, thin, dashed, dot-dashed lines correspond to 

the values for 4 ^/2a'H = l/2;l;3/2;2;5/2, respectively. Note that, in all the cases the acceleration 
parameters result to be negative. 

where y is a dimensionless function of the scalar field defined by y = aH^ (1 -|- /30)^, 7 
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is a dimensionless constant given by 7 = I I a {H()(3) and A^e corresponds to A^, 



that the energy density wil 



i./rT2^+^(i^)tanh-(VTa7. 

Let us now to consider the attractor behavior of this model. By taking into account Eq. 

(2.13), we get 

dHi<P)=6Hi<pjexp\^ [\{HoAA, (5.7) 

where 4>^ represents the initial value of the scalar field (j). The function g^Ho) is given by 

/ (1 — 2aHQ^'^ and it is positive for 2aHQ < 1 (this makes sure 

be positive, as we can see from Eq. (5.1)). Thus, the integrand 
within the exponential term will be negative, due to that d(j) and Hq have contrary signs 
(assuming that the perturbation 6H does not change the sign of (^)[17]. In this way, all the 
linear perturbations tend to vanish rapidly. 

In what concern to the scalar perturbations, the scalar spectral index parameter becomes 

n,-l = 2,„-4(l-^-^j.„ (5.8) 

and the running scalar spectral index results to be 



where ^„ is defined as 




(l-2a/7^)^^. (5.10) 



Analogously, from Eq. (3.15) we find for the gravitational wave spectral index, n^,, an 
expression given by 

n^ = -2e^, (5.11) 
and we find the corresponding tensor-to-scalar amplitude ratio 

r = \ {\-2ali^)e^. (5.12) 

Bearing in mind the expression that we have required for H{(j)), we obtain a relationship 
between r and given by 

= <|±i^, (5.13) 

I67 + 1 - Hs 

where the dimensionless constant 7 was defines previously. Note that we need to satisfy 
Us < 1 + I67 in order to have r > 0. Thus, from this inequality we get a constraint on the 
parameter 7 given by 7 > jqHus — 1)|- 

Figure 3 shows how changes r as a function of for two different values of the parameter 
7. These values are 7 = 8.0 x 10^^ and 7 = 16.0 x 10^'^. From this figure we see that our 
model can accommodate quite well the observational data. Note that this model allows the 
possibility of having a Harrison-Zel'dovich spectrum, i.e. rig = 1, with r 7^ as could be seen 
from this plot. 
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Figure 3. This plot shows the parameter r as a function of the scalar spectral index Ug for two 

values of the constant 7 = j a{Hol3) , as described by Eq. (5.13). Here, we have taken the 

values 7 = 8.0 x 10^"^ and 7 = 16.0 x 10^'^. Note that we could have the possibility of having a 
Harrison-Zel'dovich spectrum (rig = 1) with r ^ 0. 



In this case the system of flow equations (4.2) is reduced to the fohowing set of equations 





2 






dN 




2aH^ 


da 


6 


^2- 




In ^~ 




2aH^ 








6aH^ 


dN 




2aH^ 



1 - 2am j 



2e«, 



(5.14) 



1 



{I > 2). 



In order to solve this set of equations we need to have H = H{N). To get this, we start by 
considering Eq. (2.8), which results N = N{(j)). Then, we need to invert this latter expression 
(if possible) to obtain (p = (pi^)- Finally, with this expression we obtain H as a function of 
N, and, by introducing this function into the flow equation we proceed to solve them. 

There exist another way for getting a relationship between H and A^. Let us assume 
that we really know the Hubble rate as a function of the scale factor, i.e., we know explicitly 
H{a). Then, since by definition dN = —din a, we get that a{N) = aee^^''~^\ where 
and A^e are the values of the scale factor and the number of e-folding at the end of inflation. 
Then, by a direct substitution of a{N) on the Hubble rate H it is obtained H{N). 

As an example of the latter approach, let us consider the model in which there is a 
smooth exit from inflation, under the so-called decaying vacuum cosmology[51]. There, it was 
found that the Hubble parameter as a function of the scale factor becomes 



H{a) = 2He 



a? + ai 



(5.15) 
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where He = H(ae)- In this case it is obtained that 

H{N) = He[l- tanh (iVg - N)] . (5.16) 

Let us solve numerically the set of equation (5.14) for the first two slow-roll parameters, e^^ 
and r]^, when is constant equal to 0.2. Fig. 4 shows the numerical solutions for these two 
slow-roll parameters. From this figure we see that the e^^ remains almost constant (closed to 
zero) for a wide range of values of N = — N . But, for < 1 it increases to the value 
of one. Actually, for N = Ng, i.e. at the end of inflation, = 1. In the same range, i.e. 

< 1, the other slow-roll parameter, ry^, decreases from a maximum value (closed to the 
point ~ 1) to its final value r]^^ ~ 1.2 at the end of infiation. For this parameter, in the 
case > 1, is it observed from the figure that it decreases lineally. 




2 4 6 8 10 



Ne-N 



Figure 4. Numerical solutions for and rj^ from the set of equations (5.14) in the case in which 
^ const. 0.2 and H{N) = H,^ [I - tanh {N, - N)]. 



5.2 The brane- world model 

As was mentioned in the introduction we consider a five-dimensional brane scenario in which 
the Friedmann equation is modified to[15] 



8tt \ ( 



(5.17) 



where A represents the brane tension. 
Expression (5.17) can be written as 



F{H) = b 



Stt 



(5.18) 
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where b is defined by b 



Prom this latter equation and Eq. (2.3) we obtain for the scalar potential 



V{^) = A 



1 



A 



9\by [l + {l)H^] 



(5.19) 



In order to obtain an explicit expression for the scalar potential we need to introduce an 
explicit expression for the Hubble parameter as a function of the scalar field. In this respect, 
we borrow the expression put forward by Hawkins and Lidsey for the Hubble parameter [52]. 



Thus, we take H{cf)) 



coth(/3(/>) 
sinh(/3<)!>) 



where /3 is a constant given by /3 



2ttC 



arbitrary dimensionless constant. 

The scalar potential and the scale factor become 



A 



V{(j)) = -(6 - C7^)csch" 



2irC 



, with C an 



(5.20) 



and 



a{t) = -bC 



4 \ 
t + t 



1/C72 



(5.21) 



respectively [52]. Two comments are in order, first we demand that C to be less that \/6 in 
order that the potential is positive definite, and second, in the expression for the scale factor 

1/2 



it is chosen that to 



3™« 
47rAC'4 



in order to have a(0) = 0. For early time it is found that 



a ~ t^^'-"^ , therefore for inflation to be realizable we need < 1. 
The amount of comoving inflation becomes in this case 



Here, x 



2nC 



N{x) 



and Xp 



In 



27rC 



inflation, which corresponds to 



sinh(xe 
sinh(x 

(j)g where <f 
sech~ 



COsh{x end) 



(5.22) 



cosh(x) 

is the value of the scalar field at end of 



2ttC 



In what concern to the attractor solutions, from Eq. (2.13) we get that 



_3_ 2 



H 



(5.23) 



The quantity in the square bracket is positive definite, thus the difference in sing between 
H' and dcf) makes the exponential negative, and therefore, the exponential rapidly tends to 
zero, showing the attractor feature. 

Returning to the previously introduced expression for the Hubble parameter, we obtain 
for the various slow-roll parameters, which result to be given by 



C2 



[l + sech2(/3 ,/,)] 



1 + 



3 + cosh(2 /3 , 



(5.24) 
(5.25) 
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and 



By using these expressions we could obtain and r which b' 




(5.26) 



n 




and 



respectively. It is not hard to show that the following relation holds 



r = 3C^ + (n^ - 1). 



(5.27) 



Now, due to observational constraint on r, which presents an upper limit, r < 0.20 (95% CL) 
from WMAP+BAO+SN[53], where SN is the Constitution samples compiled in Ref. [54], 
and since = 0.963 ± 0.012 (excluding the Harrison-Zel'dovich spectrum in a value greater 
than 3(t)[55], we get that the parameter C should satisfy the upper bound < 0.079 it 0.004 
in order to be in agreement with the observational data. 

On the other hand, the consistency condition in this case becomes r = . ^ 

where n^, results to be = — [1 + sech(/3 0)]. Note that in the limit in which b — > oo 
we obtain the standard results r = —2n^. 

In what concern to the hierarchy slow-roll parameters equations we find 



Following an approach analogous to the previous subsection we solve this set numerically in 
the case in which the parameter remains constant equal to 0.2, and we use expression 
(5.16) for the dependence of the Hubble parameter as a function of the number of e-folding. 
The result is shown in Fig 5. Note that rj^ increases enormously close to the end of inflation. 
With this parameter much greater that one and reaches the value equal to one at the end 
of inflation, the slow-roll approximation becomes unsustainable at the end of inflation. 

6 Conclusion 

We have studied inflationary universe models in terms of a single scalar field. We have applied 
the exact solution approach to the modified Friedmann equations. After describing the main 
characteristics of the infiationary model in general terms, we described in some details two 
specific models. First, we have studied a model characterized by a modified Friedmann 

2 4 ( 3"ip, \ 

equation of the type —aH = — — p. Here, it was described the kinematical evolution 




(/ > 2). 



(5.28) 
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3.5 



3.0 - 




Figure 5. Numerical solutions for and 77^ from the set of equations (5.28) in the case in which 
^ = const. = 0.2. Here, was used that H{N) = He[l- tanh {N^ - N)]. 



in the case in which the Hubble parameter evolves as H{(j)) = Hq{1 + (3<j)). With this at hand, 
we could obtain the scalar potential, the corresponding number of e-folding and the attractor 
feature of the model. For some values of the parameters that enter into the scenario, we were 
able to characterize inflationary universe models. 

In what concern to the scalar and tensor perturbations we calculated the scalar and 
tensor power spectrum generated by the quantum fluctuations of the scalar and the gravita- 
tional fields. We determined the scalar and tensor spectrum indices in term of the so-called 
slow-roll parameters, e^^, r]^ and From these quantities we were able to write down 
explicit expressions for the different parameters. Moreover, the shape of the contours in the 
r — Us plane results to be in agreement with those given by the WMAP 7. In fact, we have 
found that the tensor-to-scalar ratio can adequately accommodate the currently available 
observational data for some values of the parameters. 

In the case of the brane-world model the functional form for the Hubble parameter 

was taken to be H{(j)) = ^J~^ sinh^^^j ' '^^lere /3 = With this expression we could 

determine all the kinematics and dynamics of the model. On the other hand, the current 
astrophysical data put an upper bound on the constant C, which becomes < 0.079ib0.004. 

One of the important point that we did not considered here was the reheating period. 
Since, in general terms, inflation is a period of supercooled expansion that, when inflation 
ends, the temperature of the universe needs to go up to a value such that it coincides with that 
corresponding to the temperature of the radiation epoch, and thus matching the Big Bang 
model. This issue, as far as we know, has not been studied under the exact approach. Perhaps, 
this study may give some insight on a deeper understanding of the period of reheating. We 
hope to address this point in the near future. 
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